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FOREWORD 


This  report  was  prepared  for  the  U.S.A.F.  by  the  Cornell  Aeronautical 
Laboratory,  Inc.,  Buffalo,  New  York,  as  Cornell  Aeronautical  Laboratory 
Report  No.  TB-495-F-14,  under  U.S.A.F. Contract  No,  W33~038-ac-17003, 
Supplemental  Agreement  No.  14.  The  contract  was  initiated  under  the 
research  and  development  project,  identified  by  RDO  No.  458-4140,  Flight 
Test  on  the  Effects  of  Down-wash  Lag  and  Aeroelasticity  on  Longitudinal 
Stability,  and  it  was  administered  under  the  direction  of  the  Aerodynamics 
Branch,  Aircraft  Laboratory,  Engineering  Division,  Wright  Air  Development 
Center,  with  Major  C.  B.  Westbrook,  H-S.A.F.,  acting  as  project  engineer. 
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-ABSTRACT 


Studies  have  been  made  to  determine  if  discrepancies  existing  between 
calculated  arid  measured  tail  loads,  responses,  and  derivatives,  as  obtained  for 
trie  F-SOA,  can  be  eliminated  or  explained.  These  discrepancies  were  evidenced 
as  (l)  a  tail  load  smaller  than  that  predicted,  (2)  the  measured  aerodynamic 
stiffness  being  significantly  greater  than  that  indicated  by  wind  tunnel  tests, 
and  (3)  an  unreasonably  large  value  of  stabilizer  effectiveness  being  obtained 
in  the  analysis  of  the  flight  data.  Certain  of  these  results  indicated  that 
investigations  should  be  made  to  establish  the  manner  in  which  aeroelasticity 
and  unsteady  flow  aerodynamics  influence  the  tail  loads  and  rigid  body  response 
of  the  F-SOA  airplane.  It  was  found  that  (a)  the  failure  of  measured  tail 
loads  to  agree  with  calculations  is  primarily  due  to  an  error  in  measurement 
due  to  limitations  in  strain  gage  location,  (b)  the  effect  of  aeroelasticity  is 
small  and  not  of  a  dynamic  nature,  (c)  the  lack  of  agreement  between  predicted  • 
and  measured  responses  at  low  Mach  numbers  is  due  to  prediction  being  based  on 
a  wind  tunnel  evaluation  of  tail-off  static  stability  which  was  different  from 
that  obtained  in  flight,  and  (d)  the  stabilizer  effectiveness  and  downwash 
cannot  be  accurately  evaluated  unless  the  tail-off  damping  in  pitch  is 
determined.  It  is  recommended  that  studies  be  made  to  develop  additional 
techniques  for  determining  the  tail-off  stability  derivatives  and  that  efforts  - 
be  directed  towards  analyzing  longitudinal  response  data  assuming  a  complex 
downwash. 
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INTRODUCTION 


This  report  presents  the  results  of  a  study  program  which  was  oriented 
to  explain  the  discrepancies  encountered  between  calculated  and  measured 
(longitudinal)  responses  obtained  from  flight  tests  of  the  F-80A  airplane.  The 
study  described  herein  is  therefore  a  continuation  of  an  extensive  program  which 
was  designed  to  investigate  the  dynamic  stability  and  aerodynamic  tail  loads  of 
an  F-80  at  high  subsonic  Mach  numbers.  Execution  of  the  program  was  the 
responsibility  of  the  Flight  Research  Department  of  the  Cornell  Aeronautical 
Laboratory  under  the  sponsorship  of  the  Aerodynamic  and  Structure  Branches  of 
the  Aircraft  Laboratory,  Wright  Air  Development  Center. 

The  reasons  underlying  the  development  of  this  theoretical  investi¬ 
gation  are  pertinent  to  the  reporting  of  this  study  and  are  presented  in  a 
following  section  of  the  report  entitled  "Background".  This  section  also 
includes  a  comprehensive  statement  of  the  problem  together  with  a  description 
of  the  approach  that  was  followed  in  obtaining  a  solution. 

The  contractual  objective  of  the  program  was  to  make  brief,  but  all 
inclusive,  studies  directed  towards  determining  -if  discrepancies  existing 
between  calculated  and  measured  responses  and  derivatives  can  be  eliminated. 
Accordingly,  consideration  has  been  given  to  means  of  improving  the  agreement  between 
predicted  and  measured  responses  and  to  means  of  improving  or  simplifying  the 
extraction  of  stability  derivatives  from  longitudinal  response  data.  Since 
certain  investigators  had  predicted  significant  aeroclastic  and  unsteady  flow 
effects  on  the  response  of  the  F-80A,  an  investigation  was  planned  to  determine 
to  what  extent  these  phenomena  would  account  for  the  observed  discrepancies. 

The  resulting  theoretical  study  examined  the  discrepancies  existing  between  the 
predicted  and  measured  tail  loads  and  in  addition,  required  that  calculations 
be  made  to  determine  the  effect  of  structural  flexibility  on  the  longitudinal 
responses  of  the  F-80A  airplane.  The  above  mentioned  work  is  fully  described 
in  this  report,  followed  by  recommendations  for  future  aeroelastic  research. 


-  1  - 
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BACKGROUND 


Dynamic  stability  and  control  research  by  means  of  direct  measurement 
of  Ihe  aircraft  response  is  a  field  in  which  the  Flight  Research  Department  of 
the  Cornell  Aeronautical  Laboratory  has  been  active  for  a  number  of  years*  The 
development  of  the  response  measurement  technique  in  research  programs  carried 
out  in  large s  relatively  slow  speed  aircraft  ultimately  led  to  an  experimental 
investigation  which  was  designed  to  measure  the  responses  and  stability 
derivatives  of  a  fighter  airplane  at  high  subsonic  M&ch  numbers,.  The  purpose 
of  this  program  was  twofold ,■  namely  — 

(1)  to  demonstrate  "that  the  response  of  a  fighter  airplane  to  a 
sinusoidal  forcing  function  can  be  measured  with  airborne  instrumentation 
without  compromising  measurement  accuracies  and 

(2)  to  investigate  the  effects  of  the  pertinent  variables  such  as 
lift  coefficient,,  C*G*  position,  Mach  number,  etc*  on  the  measured  responses 
and  stability  derivatives  for  comparison  with  wind-tunnel  results  and  the 
applicable  theory* 

By  the  time  the  F-80A  flight  program  was  about  to  get  under  way, 
other  organizations  and  investigators  had  become  actively  interested  in  the 
problem  of  aircraft  response  prediction  and  measurement*  One  of  the  questions 
which  had  been  raised  was;  were  the  original  simplifying  assumptions  made  in 
the  earlier  work  at  Cornell  still  justifiable  in  view  of  the  higher  dynamic 
pressures  and  Mach  numbers  which  would  be  obtained  with  the  F-80A*  These 
assumptions,  which  previous  flight  measurements  indicated  to  be  valid,  were  as 
follows? 


(1)  The  effects  of  frequency  upon  the  stability  derivatives  aie 
unimportant  in  the  range  of  frequencies  normally  considered,  i*e*  0  S  8 
radians  per  second* 

(2)  The  downwash  lag  can  be  approximated  by  the  time  required  for 
the  flow  to  travel  from  the  wing  to  the  tail0 

(3)  The  airplane  is  non-elastic,  thus  making  rigid  body  dynamics 
applicable  to  the  solution  of  the  response  prediction  problem* 

In  an  effort  to  obtain  theoretical  verification  for  the  first  tw© 
assumptions,  the  study  described  in  reference  (l)  was  undertaken  to  determine 
whether  the  contribution  of  non-uniform  flow  effects  is  significant  in  the  range 
of  frequencies  used  in  the  flight  test  procedure*  This  study,  which  was  mad* 
in  a  rigorous  and  thorough  manner  and  applied  numerically  to  an  P-80A,  indicated 
that  frequency  effects  and  higher  order  aerodynamic  derivatives  are  not  im¬ 
portant  in  the  range  of  frequencies  that  describe  the  flight  path  response  of 
the  rigid  airplane*  Use  of  the  exact  unsteady  flow  theory  did  result  in  a 
conoeption  of  downwash  lag  which  is  different  from  the  time  lag  concept  used  in 
the  simple  theory*  Responses  obtained  with  the  use  of  this  phase  lag  approach 
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have  a  slightly  reduced  resonance  peak  and  in  effect  look  as  if  the  damping  is 
somewhat  greater  than  is  predicted  by  the  simple  theoryc  It  should  be  em¬ 
phasized  that  this  difference  is  small  and  that  experimental  verification  should 
be  difficult  to  obtain,,  F-80A  responses  were  also  computed  by  Walkowicz  (see 
reference  2)  in  which  he  attempted  to  include  both  aeroelastic  and  non¬ 
stationary  flow  effectSo  These  calculations  do  not  indicate  the  effects  of 
each  phenomenon  separately  but  at  zero  frequency,  where  non-station ary  flow  is 
non-existant,  the  effect  of  aeroelasticity  is  seen  to  be  quite  large.  In  view 
of  the  agreement  which  had  been  previously  obtained  between  flight  measurements 
and  the  simple  theory,  considerable  doubts  were  had  as  to  the  validity  of  these 
calculations,,  Nevertheless,  they  would  either  be  substantiated  or  disproved 
by  the  results  of  the  F=80A  flight  test  program,, 

Upon  completion  of  the  F=80A  longitudinal  flight  tests  it  was  seen 
that  the  normal  acceleration  and  pitching  velocity  responses  were  basically  in 
agreement  with  what  is  predicted  by  rigid  body  and  steady  flow  aerodynamic  theory,. 
Certainly  there  were  discrepancies,  but  they  were  in  no  wise  comparable  to  the 
results  predicted  in  reference  (2).  On  the  other  hand,  a  sizable  discrepancy  was 
noted  between  the  predicted  tail  load  and  the  measured  horizontal  tail  load. 
Moreover,  this  difference  appeared  to  be  similar  to  the  effect  indicated  in 
reference  (1)  i.e.  a  reduction  of  the  peak  response,  which  resulted  from  treat¬ 
ing  the  downwash  lag  in  an  exact  manner,.  The  resolution  of  this  problem  - 
namely,  to  what  extent  does  aeroelasticity  and  non-stationary  flow  influence 
the  flight  path  and  tail  load  response  of  the  F-80A  and  to  what  extent  are  any 
or  all  of  the  observed  discrepancies  between  measurement  and  prediction  caused 
by  these  two  phenomena  -  was  of  basic  importance  at  this  time.  Unfortunately, 
the  test  program  and  dat&  analysis  had  encountered  a  large  number  of  delays; 
and  time  and  money  were  not  available  to  carefully  examine  the  experimental 
results  in  light  of  the  various  theoretical  predictions  which  had  been  made. 
Calculations  were  made  to  check  the  influence  of  unsteady  flow  effects  on  the 
horizontal  tail  load  of  the  F-80A  airplane,  but  the  results  of  these  calcu¬ 
lations  did  not  explain  the  apparent  discrepancy  between  the  measured  tail 
loads  and  the  theoretical  predictions. 

The  study  program,  presented  herein,  was  developed  primarily  in  an 
attempt  to  resolve  the  overall  situation  described  above.  Since  the  analysis 
of  the  F-80A  longitudinal  data  proved  to  be  rather  difficult  and  somewhat 
disappointing  in  certain  respects,  it  was  felt  that  a  definite  need  existed  for 
examining  the  state  of  the  art.  In  this  respect,  it  was  believed  that  the  study 
should  be  performed  with  a  completely  open  mind  as  to  the  causes  of  any  and  all 
disagreements  between  measurement  and  theory.  In  a  technical  investigation, 
such  as  this,  very  often  the  obvious  causes  are  overlooked  in  a  search  for  more 
sophisticated  answers.  It  was  felt  that  this  tendency  should  be  avoided  and 
that  efforts  should  be  made  to  establish  that  apparent  discrepancies  are  real 
and  valid. 
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In  accordance  with  this  philosophy,  it  is  believed  that  any  failure 
of  theoretical  predictions  to  agree  with  measurement  is  the  result  of  one  or 
more  of  the  following  three  situations* 

(a)  stability  derivatives  used  in  the  theoretical  calculations 
are  not  in  sufficiently  close  agreement  with  actual  derivatives; 

(b)  experimental  measurements  are  incorrect; 

(c)  the  assumption  -  that  second  order  aerodynamic  derivatives, 
unsteady  flow  effects,  and  structural  flexibility  are  negligible  for  the  F-80A  - 
is  not  valid. 

These  items  are  listed  in  the  order  in  which  they  should  logically  be 
eliminated  as  possible  causes  of  the  discrepancies  in  question.  In  practice, 
however,  the  process  of  elimination  is  not  very  orderly  when  the  researcher 
finds  himself  groping  in  the  darko  In  the  ensuing  discussion,  the  results  of 
this  investigation  are  presented  mainly  in  the  order  in  which  they  were  obtained, 
although  there  has  been  some  regrouping  for  the  purpose  of  presenting  a  more 
coherent  argument. 

The  tail  load  discrepancy  is  examined  first,  followed  by  a  calcu¬ 
lation  of  the  effects  of  aeroelasticity  on  the  response  of  the  F-80A  airplane. 
This  section  is  succeeded  by  a  discussion  of  the  factors  which  caused  the 
observed  disagreement  between  the  estimated  and  measured  response  and  con¬ 
clusions  are  reached  regarding  the  necessity  of  including  unsteady  flow  and 
aeroelastic  phenomena  in  the  prediction  of  the  F-80A  longitudinal  response. 
Finally,  in  the  last  and  major  part  of  the  theoretical  discussion,  efforts  are 
directed  towards  examining  the  discrepancies  existing  between  predicted  and 
measured  derivatives  and  the  reasons  for  same.  Consideration  is  also  given  to 
means  of  improving  and/or  simplifying  the  extraction  of  stability  derive tives 
from  longitudinal  response  data. 
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EXAMINATION  OF  THE  DISCREPANCY  EXISTING 


BETWEEN  PREDICTED  AND  MEASURED  TAIL  LOAD 


As  was  indicated  previously,  -the  lack  of  agreement  between  the  calcu¬ 
lated  and  flight  measured  horizontal  tail  load  was  the  particular  discrepancy 
which  was  primarily  responsible  for  questioning  whether 

(a)  there  are  significant  aeroelastic  effects  on  the  F-80A  airplane. 

(b)  the  lag  in  downwash  at  the  tail  is  significantly  airferent  from 
the  manner  in  which  it  is  treated  in  writing  the  simplified  equations  of  motion. 

This  discrepancy  can  be  seen  in  figures  1  through  17,  which  have  been 
taken  from  reference  (3)  for  the  purpose  of  supplementing  this  discussion. 

Note  that  at  those  values  of  frequency  where  peak  response  is  predicted,  the 
amplitude  of  the  measured  horizontal  tail  load  lies  well  below  the  calculated 
curve  whereas  measured  and  predicted  values  of  phase  angle  are  generally  in 
excellent  agreement.  Note  also  that  the  measured  normal  acceleration  and  pitch¬ 
ing  velocity  responses  are  often  in  disagreement  with  the  predicted  response, 
especially  at  low  frequencies.  Since  there  is  disagreement  between  the  measured 
and  predicted  motion  responses,  there  should  be  similar  disagreement  between  the 
measured  and  predicted  tail  loads.  This  means  that  it  is  not  relevant  to  com¬ 
pare  observed  tail  loads  with  predicted  loads,  but  rather  they  should  be  compared 
with  the  loads  necessary  to  produce  the  observed  responses  of  normal  acceleration 
and  pitching  velocity. 

The  horizontal  tail  loads  can,  in  fact,  be  written  as  a  function  of 
the  normal  acceleration  and  pitching  velocity  responses.  The  incremental  change 
in  tail  angle  of  attack  can  be  expressed  as 

fkk-  )  . 


(1) 


The  normal  acceleration  is 


On  substituting  (2)  in  (l),  joo  for  ,  and  dividing  (l)  by  % 

obtain  ** 


^Lod 

~T 


-  =  0-  ffi+  ^  ^  A 


j-  j )  co  j-  u>  dL  eLi 

+  J)  64  ^  20»{jJ.ur  ti  ^ 


where  the  subscript  "obs."  designates  the  measured  or  observed  value  of  the 
response,,  The  tail  load  is  then  given  by 


—  ^4- 

I 


Examination  of  equations  (3)  and  (4)  indicates  that  the  tail  load 
can  be  calculated  for  each  frequency  test  point  using  the  appropriate  experi¬ 
mental  values  for  Cu  ,  /^/jf  ,  and  ^  <>  These  confutations  do  require, 

however,  that  numerical  values  be  assumed  for  , 

and  o  (A  slightly  different  approach,  using  equation  (l)  directly,  would 

require  only  the  estimation  of  d  £  ,  JL  ,  and  . )  For  the 

confuted  results  shown  plotted  in  figures  18  through  21,  wind  tunnel  values 
were  used  for  all  parameters  except  J*/  ,  which  was  assumed  to  be  as  obtained 
in  the  analyses  performed  in  reference  (3)0  This  latter  assumption  was  made  in 
order  to  be  conservative,  in  that  a  smaller  produces  a  larger  tail  load. 

The  real  meaning,  of  w conservative",  as  used  here,  will  become  evident  later  in 
the  section  devoted  to  means  of  extracting  stability  derivatives  from  flight 
data. 

Figures  18  through  21  present  the  amplitude  and  phase  of  the  hori¬ 
zontal  tail  as  obtained  experimentally  and  as  confuted  by  means  of  equations 
(3)  and  (4)o  It  is  seen  that,  while  there  is  excellent  agreement  in  phase  angle, 
differences  in  computed  and  experimental  amplitudes  are  significant  and  increase 
appreciably  with  increase  in  Mach  number.  It  is  also  seen  that  the  percentage 
difference  between  the  computed  and  experimental  amplitudes  decreases  with 
increasing  frequency  of  oscillation,  except  at  0.3  Mach  number  where  the 
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difference  between  experiment  and  calculation  is  relatively  small* 

This  disagreement  between  the  experimental  tail  load  and  that  load 
which  is  indicated  by  the  observed  responses  is  in  no  way  similar  to  the  effeot 
which  was  shown  in  reference  (l)  to  be  due  to  an  exact  treatment  of  the  down- 
wash  lag.  In  fact  this  amplitude  difference  is  so  large  at  the  high  Hach  num¬ 
bers*  that  any  attempt  to  seek  a  rigorous  aerodynamical  explanation  appears  to 
be  somewhat  ludicrous.  In  view  of  the  excellent  agreement  in  phase  angle,  one 
is  inclined  to  conclude  that  inoorrect  measurement  of  the  tail  load  is  responsible 
for  this  unusual  situation,  wherein  incorrect  amplitude  values  are  obtained  with 
phase  angles  whioh  are  substantially  correct.  An  error  in  static  sensitivity 
would  produce  the  reduced  amplitude  measurements,  but  this  would  not  explain  the 
observed  variation  in  amplitude  discrepancy  with  Mach  number.  On  noting  that  the 
dimensional  load,  in  pounds,  increases  with  Mach  number  due  to  increased  dynamic 
pressure,  it  appears  that  one  possible  explanation  is  a  non-linear  relationship 
between  load  md  strain  gage  output.  Another  possible  oause,  though  an  uncertain 
faotor,  is  the  relieving  of  tail  load  caused  by  fuselage  bending  or  other  aero- 
elastio  effects.  This  is  amplified  further  in  the  following  seotion  where  it  is 
demonstrated  that  the  responses  of  the  F-80A  are  somewhat  reduoed  by  aero- 
elasticity.  Since  aeroelastic  effects  are  a  funotion  of  dynamio  pressure,  this 
phenomenon  may  possibly  be  a  faotor  here.  Admittedly  this  is  a  speculation,  but 
in  view  of  the  fact  that  the  phase  angles  agree  at  0.75  Maoh  number  while  a  40 
percent  difference  in  amplitudes  exists  at  low  frequencies,  it  appears  that  a 
measurement  error,  in  possible  conjunction  with  aeroelastic  effects,  is  the  only 
means  of  explaining  this  result. 

The  variation  in  percent  amplitude  difference  with  frequency  oan  be 
explained  as  a  function  of  change  in  spanwise  load  distribution  with  frequenoy. 

For  example ,  reference  3  points  out  that  incremental  loads  inboard  of  station  10 
on  the  horizontal  stabilizer  were  not  measured  by  the  strain  gage  installation 
on  the  F-80A  airplane.  The  area  inboard  of  station  10  1b  16.9  peroent  of  the 
total  horlsontal  tail  area.  If  the  assumption  is  made  that  the  spanwise  load 
distribution  contributed  by  the  elevator  is  measured  in  toto  and  that  the 
measured  loading,  contributed  by  4  ,  c/.  p  and  dc  ,  is  in  error  by  sixteen 
peroent,  we  obtain  the  result  shown  in  figure  22.  Note  that  little  change  in 
phase  occurs  but  that  a  reduction  in  amplitude  does  result  in  a  manner  similar 
to  that  obtained  in  the  F-80A  flight  tests.  It  then  beoomes  feasible  to  explain 
the  variation  in  amplitude  discrepancy  with  frequenoy  by  a  measurement  error 
similar  to  that  assumed  above,  Binoe  at  the  lower  frequencies  the  greater 
portion  of  the  tail  load  is  contributed  by  &  ,  ,  and  <&.  while  at  the 

higher  frequencies  the  tail  load  is  due  mainly  to  elevator  deflection. 
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AERCBLASTIC  EFFECTS 


OH  THE  RESPCHSE  OF  THE  F-80A  AIRPLAHE 


While  the  F-80A  flight  test  data  fail  to  corroborate  the  large  aero- 
elastio  effeots  reported  in  referenee  (2),  there  is  suffioient  Justification 
for  a  quantitative  aeroelastio  study  to  be  made  as  a  part  of  this  overall 
investigation.  Such  a  study  would  indicate  the  actual  magnitude  of  the  effeots 
to  be  expeoted  for  a  relatively  stiff  airplane.  In  addition,  these  data  would 
serve  to  substantiate  the  observed  responses  in  demonstrating  that  the  F-80A 
airplane  is  not  suitable  as  an  aeroelastio  researoh  test  vehiole. 

In  performing  these  aeroelastio  oaloulations,  the  assumption  was 
made  (as  was  done  in  reference  2)  that  struotural  defleotions  of  the  wing  would 
be  negligible  and  that  only  fuselage  bending,  stabiliser,  and  elevator  twist 
need  be  considered.  Stabiliser  and  elevator  twist  were  not  properly  included 
as  additional  degrees  of  freedom  in  the  rigid  body  equations  of  motion,  as  was 
fuselage  bending,  because  of  the  complexities  involved.  Rather  the  computed 
statlo  reduction  in  elevator  effectiveness  (stabiliser  twist  was  found  to  be 
negligible)  was  as  Burned  to  be  valid  over  the  entire  frequency  range  sinoe  the 
dynamios  of  the  elevator  structure  are  far  removed  from  the  frequenoy  range  of 
interest.  This  assumption  of  constant  effectiveness  is  substantially  correot, 
if  the  major  part  of  the  loading  that  onuses  the  elevator  to  twist  is  caused  by 
elevator  deflection,  as  is  presumed  to  be  the  case.  The  method  used  in  computing 
the  reduced  sero  frequenoy  value  of  elevator  effectiveness  is  that  desoribed  in 
reference  (4).  For  the  configuration  selected  for  oaloulation  purposes  (Haoh 
numbers  .6,  Altitude!  10,000  feet),  a  seven  peroent  static  reduction  in  elevator 
effectiveness  was  found  for  the  F-80A  airplane. 

In  introducing  the  additional  degree  of  freedom  into  the  equations  of 
motion  for  the  purpose  of  representing  a  flexible  fuselage,  it  was  assumed  lhat 
the  fuselage  bending  was  in  phase  with  the  aerodynaado  and  inertia  foroes. 

Stated  differently,  the  dynamics  of  the  fuselage  structure  were  negleoted  sinoe 
\ he  fuselage  bending  modal  frequenoy  is  presumably  much  higher  than  the  rigid 
body  frequenoy  range  under  consideration.  This  simplifying  assumption  permits 
the  dlreot  formulation  of  the  fuselage  bending  equation* 


where  Jf  is  the  fuselage  bending  angle  measured  as  the  change  in  angle  of 
attaek  of  the  horizontal  stabiliser. 
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is  the  coefficient  of  fuselage  bending  due  to  horizontal' tail  load* 
is  the  coefficient  of  fuselage  bending  due  to  normal  acceleration  of  the 


o.g. 


is  the  coefficient  of  fuselage  bending  due  to  pitching  acceleration* 


The  coefficient  of  fuselage  bending  due  to  horizontal  tail  load  was  cbeoked 
experimentally  and  found  to  agree  exceedingly  well  with  the  value  computed  in 
reference  (2).  In  view  of  -this  excellent  agreement,  .  a  value  for  jj)  jT 

was  also  taken  directly  from  this  source*  By  working  backwards  from  the  bend¬ 
ing  moment  curve  due  to  fuselage  dead  weight,  the  fuselage  weight  distribution 
curve  was  found .  It  was  then  possible  to  compute  the  fuselage  bending  moment 
due  to  angular  acceleration  about  the  o.g.  and  subsequently  determine  the 
ooeffioient  of  fuselage  bending  due  to  pitching  acceleration.  The  results 
obtained  were  as  follows  t 

*=/.</ x/o-‘ 


|£-  =  -/«?/  X/D-’ 

j2£.  =  _  4.2 S X /O'*  sect 
?  6 


and  JT 


Equation  (5)  can  be  written  in  terms  of  the  variables  &£  , 
by  noting  that 


DO  *  %  , 


i-^i) 


(6) 


and 


(7) 
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On  substituting  equations  (6)  and  (7)  into  (5),  we  obtain 


=  (a.  D  fa)  CxL  -i  (c  O  +  £ )p&  -j.  /  $  (8) 

where 
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The  longitudinal  equations  of  motion,  with  the  longitudinal  (x)  degree  of  free¬ 
dom  neglected  then  becomes 


c< 

f 

_ l 

?„)  (teD'rn^)  -v^/77^ 

pa 

= 

-(a~Djjb)  -(cDtc)  / 

r  _ 

/i 

where 

j.  )  — 

and 


We  find  therefore,  under  the  assumptionB  expressed  in  equation  (5), 
that  the  order  of  this  system  of  equations  is  still  twoc  This  is  as  was  to  be 
expeoted  since  the  structural  transfer  function  of  the  fuselage  has  been  assumed 
to  be  unity,  which  assumption  is  believed  to  be  in  accordance  with  the  facts . 

Use  of  the  above  equations  of  motion  permit  the  comparison  of  responses  obtained 
for  both  a  rigid  airplane  and  one  possessing  a  flexible  fuselage „  Figures  23 
through  25  present  the  amplitude  and  phase  of  the  pitching  velocity,  normal 
acceleration  and  horizontal  tail  load  responses  of  the  F-80A  for  both  an  elastic 
and  rigid  body  computationo  Note  that  little  or  no  difference  in  phase  angle  is 
obtained  in  the  responses  of  the  elastic  and  non-elastic  airplane,,  The  reduction 
in  response  amplitude  appears  to  be  similar  to  what  would  result  from  a  straight 
reduction  in  elevator  effectiveness.  This  is  not  actually  the  case,  however, 
since  the  fuselage  deflection  response  varies  with  frequency  (see  figure  26) 
and  therefore  the  percentage  reduction  in  response  amplitudes  does  likewise,  A 
plot  of  the  peroent  change  in  amplitude  caused  by  the  combined  effect  of  re¬ 
duced  elevator  effectiveness  and  fuselage  bending  is  presented  in  figure  27 
for  the  three  responses  shown  in  figures  23  through  25.  Since  at  steady  state 
(zero  frequency)  a  7  percent  reduction  is  caused  by  elevator  twist,  the  re¬ 
maining  9  percent  reduction  in  the  pitching  velocity  and  normal  acceleration 
responses  is  oaused  by  fuselage  flexibility. 
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Experience  indicates  that  it  would  be  extremely  difficult  to  detect 
the  above  calculated  aeroelastic  effects  in  the  results  of  an  actual  experi¬ 
mental  flight  program.  The  nature  of  the  computed  effeot  is  suoh  that  for  the 
F-80  airplane,  ae roe las ti city  aots  more  like  a  scalar  factor  instead  of  a 
vector  quantity.  Moreover,  the  average  value  of  this  soalar  factor  falls  into 
the  region  of  prevailing  prediction  and  measurement  errors,  A  careful  exam¬ 
ination  of  the  response  plots  presented  in  figures  1  through  17  fails  to  produoe 
any  positive  verification  of  these  computed  aeroelastic  effects. 
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C0ISIDKRATI08  OF  MBiHS  OF  IMPROYISG  THK  AQRHMEHT 


BKTTISEN  PREDICTED  AMD  MEASURED  LOHGITPDIHAL  RESPQ1SIS 


Ap  examination  of  the  F-80A  response  plots  presented  in  figures  1 
through  17  indicates  that  the  measured  responses  are  representative  of  a 
second  order  dynamic  system.  On  the  basis  of  this  visual  comparison  between 
calculation  and  experiment,  it  cannot  be  determined  whether  this  system 
possesses  oonstant  derivatives  or  parameters  which  vary  with  the  frequency  of 
excitation.  It  is  pointed  out  in  reference  (1)  that  the  differences  in  response 
calculations  obtained  with  the  use  of  the  simple  and  exact,  complex  equations 
of  motion  are  so  small,  that  the  inaccuracies  of  using  the  time-lag  oonoept  of 
downwash  are  no  greater  than  those  of  the  original  theoretical  assumptions  or 
the  experimental  data.  Thus,  comparison  of  response  data  is  not  a  valid  means 
for  evaluating  the  effects  of  oertain  phenomena  when  these  effects  are  as  small 
or  smaller  than  the  errors  which  are  introduced  by  faulty  measurement  or  faulty 
prediction.  This  resulting  inability  to  detect  these  effects  in  this  manner 
does  not  say  that  they  are  not  present,  however.  On  the  other  hand,  if  one 
attempts  to  account  for  these  phenomena  in  the  analysis,  which  is  performed  to 
determine  the  stability  derivatives,  one  finds  that  one  must  deal  with  a  com¬ 
plicated  system  whose  analysis  is  extremely  difficult  or  perhaps  impossible. 

In  the  following  seotion  of  this  report  consideration  is  given  to  the 
analysis  part  of  the  above-mentioned  problem.  As  far  as  improved  agreement 
between  response  prediction  and  measurement  is  oonoemed,  it  appears  that  major 
attention  should  be  first  directed  towards  increasing  the  accuracy  with  which 
the  longitudinal  short  period  natural  frequency  and  damping  are  originally  esti¬ 
mated.  Certain  of  the  F-80A  response  plots  (particularly  those  for  the  23  per- 
oent  and  27  peroent  o.g.  positions)  indieate  that  the  aerodynamic  stiffness  of 
the  airplane  is  considerably  different  from  what  was  originally  predicted.  If 
there  are  differences  between  the  flight  measured  and  predicted  values  of  damp¬ 
ing  for  the  airplane,  the  response  plots  indicate  that  these  differences  are 
not  as  large  or  as  significant  as  the  discrepancies  encountered  between  the 
measured  and  predicted  values  of  stiffness  or  natural  frequency. 

The  above  observations  point  out  the  necessity  for  accurately  determin¬ 
ing  the  static  stability  or  neutral  point  of  the  airplane  when  attempting  to 
prediot  its  longitudinal  response.  In  the  oase  of  the  F-80A,  these  predictions 
were  based  on  wind  tunnel  data  whioh  were  obtained  from  reference  (5).  The 
large  variation  in  neutral  point  or  /<?£C 1  shown  by  the  dashed  line  in 

figure  28  is  the  result  of  the  non-linearity  found  to  exist  in  the  pitching 
moment  data  obtained  from  the  wind  tunnel  tests.  This  non-linearity  necessi¬ 
tated  that  the  slope  of  the  pitching  moment  curve  be  determined  for  the 
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particular  lift  coefficient  and  Mach  number  configuration  at  which  it  was  in- 
tended  to  test  tiis  airplane  in  flight.  It  is  seen  that  any  deviation  in 
airplane  weight  and  resulting  lift  coefficient  in  flight  from  that  selected 
for  calpulation  purposes  will  automatically  cause  a  small  discrepancy  in 

o  The  non-linearity  of  these  pitching  moment  curves,  therefore, 
complicates  the  task  of  accurately  predicting  the  flight  neutral  point. 
Examination  of  the  wind  tunnel  da^a  shows  that  this  non-linearity  is  due  en¬ 
tirely  to  the  pitching  moment  characteristics  of  the  model  minus  the  horizontal 
tail.  Since  the  contribution  of  the  horizontal  tail  to  the  static  stability,  in 
conjunction  with  wind  tunnel  measured  values  of  stabilizer  effectiveness  leads 
to  wind  tunnel  values  of  downwash  which  are  in  excellent  agreement  with  downwash 
theory,  there  is  reason  to  believe  that  the  tail-off  static  stability  data  is 
primarily  responsible  for  the  observed  disagreement  in  predicted  and  measured 
stiffness. 


This  hypothesis  is  substantiated  by  the  results  obtained  in  static 
flight  tests  made  with  the  F-80A.  Reference  (6)  shows  that  the  neutral  point 
of  the  F-80A  is  significantly  farther  aft  (at  high  lift  coefficients)  than  what 
was  indicated  by  the  above-mentioned  wind  tunnel  tests.  This  result  intimates 
that  the  extreme  non-linearity  noted  in  the  wind  tunnel  test  data  is  in  error 
and  that  the  airplane  is  considerably  more  stable  at  high  angles  of  attaek  than 
is  shown  by  the  dashed  line  in  figure  28.  Since  it  has  been  noted  that  knowledge 
of  the  static  stability  is  critical  for  an  accurate  prediction  of  the  longitudinal 
frequency  response,  it  would  seem  advisable  that  static  flight  tests  be  made, 
whenever  possible,  to  check  wind  tunnel  results.  These  static  tests  would  per¬ 
mit  -the  accurate  evaluation  of  the  elevator  power,  ,  in  addition  to  the 
location  of  the  neutral  point.  Admittedly,  flight  tests  are  not  a  convenient 
means  of  investigating  -the  effects  of  Mach  number,  but  in  this  particular 
instance,  it  is  the  variation  of  static  stability  with  angle  of  attack  that  is 
question  able. 


Note  that  static  flight  tests  should  be  performed  in  the  glide  config¬ 
uration,  in  order  to  prevent  the  direct  thrust  effects  from  contributing  to  the 
static  stability.  In  constant  speed  flight,  which  is  the  case  for  the  short 
period  mode,  the  jet  thrust  does  not  contribute  to  the  airplane's  stability.  It 
does,  however,  contribute  to  the  static  stability  when  the  speed  is  allowed  to 
vary.  This  is  demonstrated  quite  readily  by  writing  the  moment  coefficient  due 
to  the  thrust  acting  at  some  vertical  arm  to  the  airplane's  center  of  gravity. 


Cv, 


f-  S  c 


where 


is  the  distance  from  the  c.g.  to  -the 


thrust  line. 


'S'1 
*  ^  / 
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If  the  airplane  is  considered  in  unaocelerated  flight  (i,e0  C u  varies  with 
speed)  the  dynamic  pressure ;  ;  will  be  the  following  function  of  lift  coef¬ 

ficient! 


Equation  (9)  becomes 

r 


\A/  c 


(10) 


For  a  constant  throttle  setting,  the  thrust,  7^  ,  is  nearly  independent  of 
airplane  speed  and  can  be  considered  a  constant.  The  stability  contribution 
oan  be  obtained  by  differentiating  (10)  with  respeot  to  Cu  . 


Although  the  direct  thrust  effect  is  not  a  factor,  a  jet  power  plant 
does  contribute  to  the  dynamic  longitudinal  stability  of  an  airplane.  These 
contributions  include  (a)  the  direot  normal  force  effects  at  the  air  duct  inlet, 
(b)  the  effect  of  the  induced  flow  at  the  tail  duo  to  the  inflow  to  the  jet 
blast,  and  (c)  an  increment  in  as  given  by  Braun  in  reference  (7),  The 

geometrical  configuration  of  the  F-80A  airplane  is  such  that  the  first  two 
effeots  are  found  to  be  negligible.  The  damping  effect  of  the  jet  stream  is 
also  quite  small  but  it  oan  be  a  significant  portion  of  the  damping  due  to  the 
-fuselage  and  the  wing.  Since  the  dimensional  damping  moment,  ^ 

is  approximately  constant  with  airspeed,  the  non-dimensional  derivative, , 
is  found  to  vary  inversely  with  airspeed  at  a  given  altitude.  For  an  altitude 
of  20,000  feet,  11  *  0,3,  it  was  oomputed  that 
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For  the  same  conditions  (based  on  wind  tunnel  data)*, 

These  numerical  results  are  presented  for  the  purpose  of  indicating  the  relative 
magnitude  of  the  damping  caused  by  the  jet  stream.  In  order  to  predict  the 
damping  of  the  short  period  as  accurately  as  possible *,  the  derivative *>  , 

should  be  carefully  estimated  by  including  all  sources  of  damping  such  as  the 
jet  damping  mentioned  above , 

In  this  discussion  of  means  of  improving  the  agreement  between  pre- 
dieted  and  measured  responses*,  it  should  be  noted  that  the  damping  in  pitch  of 
the  F-80A  was  assumed*,  for  calculation  purposes*,  to  be  caused  only  by  the  hori¬ 
zontal  tail.  In  other  words*,  the  tail  off  damping  in-  pitch  was  assumed  to  be 
zero.  While  this  rough  approach  to  the  estimation  of  the  derivatives  may  be 
considered  oversimplified  from  the  aerodynamicists’  point  of  view*  it  can  be 
justified  engineering-wise  since  the  response  predictions  were  desired  pri¬ 
marily  for  the  purpose  of  instrumentation  selection  and  to  provide  a  check  on 
flight  measurements. 

It  is  believed  that  the  above  discussion  sufficiently  emphasizes  the 
point  which  was  made  earlier!  namely  -  the  errors  introduced  by  the  estimation 
of  tiie  first  order  aerodynamic  terms  are  large  enough  to  completely  mask  the 
effects  of  the  second  order  terms,  the  unsteady  flow  effects*,  and  the  possible 
aeroelastic  effects  for  the  F-80A  airplane.  From  the  response  prediction 
standpoint*,  consideration  of  these  latter  items  for  the  F-80A  airplane  appears 
to  be  purely  academic. 
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EXAMINATION  OF  AND  MEANS  FOR  IMPROVING  AHP/OR  SIMPLIFYING  THE  EXTRACTION 


OF  STABILITY  DERIVATIVES  FROM  LONGITUDINAL  RESPONSE  DATA 


1,  General  Considerations. 

Before  proceeding  to  analyze  response  data  obtained  in  flight  it  is 
necessary  to  assume  equations  of  motion  which  adequately  describe  the  dynamics 
of  the  system.  Past  experience  has  shown  that  linear,  second  order  differential 
equations  mathematically  represent  the  short  period  longitudinal  mode  of  a 
rigid  airplane  with  excellent  accuracy  even  when  the  limitation  of  infinitesimal 
displacements  is  violated.  In  this  case,  the  assumption  of  linearity  and  the 
restriction  to  first  order  aerodynamic  and  inertia  forces  not  only  eases  the 
prediction  problem  but  simplifies  greatly  the  problem  of  analysis,  which  has  as 
its  objective  the  extraction  of  stability  derivatives  from  experimentally 
obtained  flight  data.  Often  there  is  some  question  whether  these  linear 
equations,  containing  a  limited  number  of  terms,  are  valid  for  certain  aircraft 
configurations  and  flight  conditions.  This  question  is  prompted  by  any  analy¬ 
sis  in  which  the  experimental  scatter  is  found  to  be  large  or  the  derivatives 
appear  to  vary  with  the  frequency  of  oscillation.  Disregarding  for  the  moment 
the  possibility  of  measurement  error,  it  becomes  incumbent  upon  the  analyst  to 
consider  other  phenomena  which  conceivably  might  influence  the  aircraft  response. 

At  this  point,  considerable  engineering  judgement  must  be  exercised  in  select¬ 
ing  the  additional  complexities  to  be  introduced  into  the  equations  of  motion. 

If  these  additional  complexities  cause  the  problem  to  become  non-linear,  it  may 
not  be  possible  to  perform  an  analysis  at  all.  Finally,  if  the  effects  of 
these  additional  phenomena  are  small  (as  they  certainly  are  for  the  F-80A)  it 
will  be  very  difficult  to  accurately  extract  the  pertinent  parameters  describ¬ 
ing  these  effects. 

Since  increasing  the  rigor  with  which  the  dynamical  system  is  represented 
simultaneously  increases  the  complexity  of  the  analysis,  it  appears  that  it  would 
be  wise  to  use  the  simplest  approach  in  the  analysis  unless  there  is  substantial 
evidence  that  unsteady  flow  or  aeroelastic  effects,  for  example,  are  present  to 
a  significant  degree.  If  this  evidence  or  knowledge  is  present,  we  soon  discover 
that  sizable  advancements  must  be  made  in  the  science  of  solving  the  reverse 
dynamics  problem  in  order  that  these  more  complicated  systems  may  be  successfully 
analyzed.  In  this  section  of  the  report  some  consideration  will  be  given  to  a 
possible  means  of  performing  an  analysis  of  longitudinal  response  data  when  the 
equations  of  motion  are  written  in  terms  of  the  complex  downwash  rather  than 
the  time  lag  concept.  The  method  will  be  indicated,  rather  than  applied,  for 
reasons  to  be  given  later.  Before  examining  the  results  obtained  in  the  original 
F-80A  analysis,  consideration  is  given  below  to  the  peculiar  problems  encountered 
in  the  analysis  of  longitudinal  response  data. 
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2.  Analysis  Methods. 

In  this  section,  discussion  Hill  be  confined  to  techniques  of  analys¬ 
ing  the  short  period  longitudinal  mode  only,  si  nee  the  experimental  flight  data 
cover  a  frequency  range  of  approximately  1  rad./seo.  to  8  rad./sec.  It  will  be 
assumed  that  the  motion  of  the  airplane  is  represented  by  linear  differential 
equations. 

a.  Operational  Analysis  Using  Tail  Load  Data  to  Determine  Basio  Derivative* 

Previous  investigators  have  shown  that  extracting  the  ouustants,  found 
in  the  operational  (longitudinal)  equations,  from  response  data  obtained  in 
flight  is  a  relatively  straightforward  task.  These  constants  are  the  transfer 
function  coefficients,  such  as  ba  «nd  fco  -  the  damping  and  stiffness  para¬ 
meters  obtained  from  the  characteristic  or  determinentftl  equations  of  the 
system.  Stability  analysts  have  also  determined  that,  in  the  ease  of  a  heavily 
d craped  short  period  mode,  the  transfer  funotion  or  operational  constants  cannot 
readily  be  obtained  from  transient  response  time  histories  of  the  system.  They 
can,  however,  be  extracted  from  frequency  response  data  (obtained  either  from 
steady  state  oscillations  or  by  harmonic  analysis  of  transient  responses)  by 
means  of  various  curve  fitting  techniques,  either  numerical  or  graphical.  The 
graphical  technique  has  evolved  from  the  methods  of  the  servomechanism  engineer 
and  proves  to  be  practical  when  the  system  is  only  second  order.  Hote  that 
these  transfer  function  parameters  are  combinations  of  the  actual  stability 
derivatives  and  as  suoh  do  not  indicate  the  basio  aerodynamic  characteristics 
of  the  airplane.  In  reference  (3)  an  analysis  procedure  was  developed  whereby 
the  tail  load  data  could  be  used  to  break  down  the  quantities,  jb0  and  , 

into  their  component  derivatives. 

On  neglecting  speed  changes  and  non- stationary  terms,  it  was  assumed 
in  the  original  F-80A  analysis  that  the  pertinent  derivatives  were  as  followst 


Since 
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where  the  subscript  *i  "  denotes  the  tail-off  derivative,  it  is  evident  that 
the  basic  quantities  to  be  evaluated  are  » 

and  ^  £  o  The  normal  acceleration  and  pitching  velocity  responses  were  put 

into  the  following  operational  form  for  convenience  of  analysiso 


ft  A 

/  4?  &o  C.  O  -3 


where 


Using  the  experimental  values  of  ,  £1  »  T" »  and and  the  measured 
amplitudes  and  phases  of  the  acceleration  and  pitching  responses,  reference  (3) 
presents  an  iterative  procedure  whereby  the  values  of  /C0  are 

evaluated  for  each  frequency  test  pointo  The  elevator  effectiveness,  /7?j  ,  is 
determined  from  the  slope  of  a  straight  line  defined  by  points  calculated  at 
all  test  frequencies* 

To  evaluate  the  remaining  unknowns,  this  procedure  was  repeated  using 
similar  expressions  based  on  tail-off  equations  of  motion*  The  derivation  of 
these  equations  is  given  below  and  it  is  seen  that  the  horizontal  tail  load 
must  be  measured  in  order  to  carry  out  this  analysis*  On  writing  the  longi¬ 
tudinal  equations  of  motion  in  terms  of  the  aerodynamic  forces  contributed  by 
the  horizontal  tail  and  the  airplane  minus  the  tail,  we  have 

&  *  +h'  06  ~ 

a?  +  mh  do  -  m^pe du  /  cg  0*0 
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where  the  subscript  "t"  denotes  the  horizontal  tail  contribution. 

Since  rj 

and 

we  obtain  on  neglecting  ,  the  following  two  equations  of  motion  in  which 

the  horizontal  tail  load  is  now  the  forcing  function  instead  of  the  elevator. 


fa,  «-  +  09- Doc  ~~jk/c  fa  (13) 


ol  +  m*,  do  -  tao!e  =  \/fa  (i4) 


The  normal  acceleration  and  pitching  velocity  responses  to  sinu^ 
soidal  tail  loads  can  be  written,  similarly  to  equations  (ll)  and  (12)  „  as 
follows*  ' 
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where 


_ 

^7-  2 to? 


As  indicated  previously,  use  of  the  above  equations,  together  with  the  measured 
responses  and  the  quantities  ia  ,  ,  7"  ,  and  ^ ,  permits  the  evaluation  of 

s  g  *nd  /f0,  »  for  each  frequency  test  point.  The  tail-off  deriva¬ 

tives,  mJ  and  ft)*/  t  are  then  found  from  definitions  of  the  tail-off  damping 
and  stiffness  parameters.  Finally,  the  last  two  unknown  parameters, 
and  ^ t  ,  are  obtained  from  a  simultaneous  solution  of  the  following  two 

oTZ. 

equationsi 


Values  of  ths  derivatives  as  obtained  by  this  method  were  plotted  against  the 
non-dimensional  frequency,  ,  and  averaged  arithmetically.  It  was  believed 

that  this  point  by  point  analysis  adequately  accounted  for  the  variations  in 
aircraft  weight  and  velocity  which  occurred  during  the  test  runs. 


b.  Equation  of  Motion  Analysis 

In  order  to  extract  stability  derivatives  directly  from  response  data 
obtained  in  flight,  it  becomes  necessary  to  curve-fit  the  data  on  the  basis  of 
the  original  equations  of  motion.  Unfortunately,  the  longitudinal  equations 
are  such  that  a  linear  dependency  exists  between  the  variables  of  motion  which, 
in  the  reverse  dynamics  problem,  are  the  coefficients  of  the  unknown  derivatives 
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Mathematically  this  means  that  a  unique  set  of  values  for  the  derivatives  cannot 
he  obtained  since  the  determinant  of  the  coefficients  (in  this  case,  the  variables 
of  motion)  proves  to  be  sero.  ThiB  is  readily  demonstrated  by  writing  the 
equations  of  motion  in  a  form  adaptable  to  their  inverse  solution  for  the  un¬ 
known  derivatives,  visj 


r 


pe 


7-  = 
% 

>. 

now 


/71* 


The  subscript  *  varies  from  1  to  /?  where  /?  is  equal  to  the  total  number  of 
unknowns.  (In. this  case,  four.)  It  is  seen  that,  in  the  inverse  solution,  the 
coefficients  of  the  unknowns  are  the  four  variables  of  motion  ,  o  ,  £>&  „ 

and  QoL  „  Since 


we  find  that 

oL  +  (^j)oa  +  o 


which  expression  illustrates  the  linear  dependency  existing  between  the  four 
quantities  (  ,  f  ,  P*L  ,  D&  )»  It  can  then  be  shown  that 
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In. reference  (8),  Schumacher  suggests  that  the  linear  dependency 
be  removed  from  the  longitudinal  equations  of  motion  by  making  the  following 
substitutions  for  /7?.  and  0 


/77* 


fpl  cl 


Note  that  this  is  an  approximation  which  assumes  that  the  tail-off  damping  in 
pitch  is  zero,  since  accurately  expressed  within  the  limitations  of  the  time 
lag  concept  of  downwash 


iv 


__  aSs 


With  the  approximate  expression  for  HI*,  ,  the  longitudinal  equations  of  motion 
(in  the  form  suitable  for  their  inverse  solution)  are  reduced  to  three  unknowns 
and  three  linearly  independent  coefficients,  viz? 


-  23  - 
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If  the  tail~off  damping  in  pitch  is  not  assumed  to  be  aero*  the  above  equations 
can  be  written  more  accurately  as 

^  S  pe  (pe  +  g%D^  ^  ny 

.  ‘  P  L  (21) 


V1"'\ 

(pe  i^poi)] 

^ms 


ie  D’e  <22> 


In  this  case  there  are  four  linearly  independent  coefficients*  namelys 


U.  ^  ,  P&,  aud  P  9  +  Pot 
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Note  that  the  subscript  "  i  *  indicates  that  one  substitutes  values  for  the 
coefficients  obtained  at  different  time  intervale,  if  he  is  dealing  with 
transient  data,  or  values  obtained  at  discrete  intervals  of  frequency,  if  he  is 
working  with  frequency  response  data.  In  the  latter  case,  eaoh  ooefficient  is 
a  vector  or  complex  quantity  and  must  be  treated  as  suoh  when  normalising  the 
equations  written  for  a  large  number  of  frequency  points  to  yield  the  same 
number  of  equations  as  there  are  unknowns.  In  reference  (8),  Schumacher 
chooses  the  method  of  least  squares  to  average  the  data.  This  normalisation 
procedure  is  indicated  in  matrix  notation,  as  follows* 


P  ■  total  number  of  test  points 


Solutions' of  equations  (19)  through  (22)  by  the  above -indicated  procedure 
yields  values  of  9  derivatives  which  are  associated  with  the  artificially 
created  motion  variable,  pg  /  ^  p ^  •  These  derivatives  will  be 

slightly  in  error  when  the  assumption  is  made  that  the  tail-off  contribution  to 
the  pitching  velocity  derivatives  is  aero.  The  accuracy  of  these  derivatives 
is  also  a  funotion  of  the  assumed  values  of  ,  but  Schumacher  feels  that 

the  effeot  of  inaccuracies  in  the  assumed  values  of  6^*  should  be  minor.  Both 

the  analysis  method  presented  in  reference  (8)  and  the  modification  suggested 
above  have  been  applied  to  some  of  the  F»80i.  flight  data  in  an  effort  to  explain 
some  of  the  discrepancies  which  are  found  to  exist  between  the  estimated  and 
measured  derivatives.  The  significant  results  are  presented  in  tte  following 
section. 
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3.  Examination  of  Discrepancies  Between  Measured  and  Predicted 
Derivatives  -  F~80A  Airplane. 

Figures  28  and  28  present  a  comparison  between  the  estimated  deriva- 
tives  (based  on  wind  tunnel  data)  and  the  measured  derivatives  as  determined 
from  the  flight  test  data  by  the  analysis  method  which  makes  use  of  tail  load 
measurements,,  Figure  30  shows  estimated  and  flight  test  determined  values  of 
the  damping  and  stiffness  parameters.,  A  and  &  .  These  curves  represent 

only  one  altitude  and  C.G.  configuration  of  the  total  flight  configurations 
represented  by  figures  1  through  17.  They  are,  however,  more  than  sufficient 
to  indicate  that  significant  differences  exist  between  prediction  and  flight 
test  re suits o 

Examination  of  figure  30  shows  that  the  flight  measured  values  of 
the  stiffness  and  damping  of  the  airplane  are  greater  than  that  predicted*  with 
the  measured  values  of  /fo  being  drastically  larger,  especially  at  Q.„3  Mach 
number,,  This  discrepancy  in  stiffness  is  likewise  reflected  in  figure  28, 
where  the  static  stability*  ^ s^lown  plotted  versus  Mach  number® 

The  fact  that  the  airplane  is  more  stable  than  what  was  indicated  in  the  wind 
tunnel  tests  is  not  the  major  cause  for  concern.  HThat  is  disturbing,  however, 
is  that  figures  28  and  29  indicate  that  this  increase  in  stability  is  caused 
primarily  by  an  increase  in  stabiliser  effectiveness.  This  increase  is  of  the 
order  of  25  to  35  percent  greater  than  those  values  of  yielded  by  the 

*•  tr 

wind  tunnel  tests  reported  in  reference  (5).  Investigation  showed  that  the 
value  of  the  tail  lift  curve  slope,  corresponding  to  these  values  of  stabilizer 
effectiveness,  would  be  aerodynamioally  impossible  to  achieve  even  if  the  tail 
efficiency  were  assumed  to  be  one  hundred  percent. 

Further  examination  of  figure  29  indicates  that  the  measured  elevator 
effectiveness,  ,  is  correspondingly  higher  than  the  value  obtained  in  the 

wind  tunnel.  The  lift  curve  slope  obtained  from  flight  data  is  approximately 
7.5  percent  less  than  the  estimated  value,  while  the  slope  of  the  downwash 
appears  to  be  unreasonably  high.  Since  /T'V  and  were  computed  from  the 
flight  measured  values  of  .  and  ct £  ,  these  derivatives  are  likewise 

*<•  W* L 

larger  than  the  original  prediction. 

In  order  to  check  these  discrepancies,  efforts  were  first  directed 
towards  an  examination  of  the  tail-on  analysis,  which  as  indicated  above  yields 
values  of  ha  ,  ,  and  ^ .  The  possible  effects  of  measurement 

errors  as  they  would  influence  results  of  this  analysis  were  considered  in  this 
phase  of  the  investigation.  These  studies  indicated  that  errors  in  static 
sensitivity  (for  a  particular  recording  channel)  would  affect  only  the  evaluation 
of  the  derivatives  /79<  and  .  Examination  of  equation  (ll)  reveals  that 

the  transfer  function  constants,  6<>  and  ^  ,  are  a  function  only  of  the 

phase  angle,  .  Further  studies  indicated  that,  whereas  reference  (3)  proposes 
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an  iterative  procedure 
the  angle,  p  ,  can  be 


for  determining  after  first  assuming  a  value  of 
expressed  empirically  with  excellent  accuracy  as 


+■  ^,0  4-  70’ 


(23) 


This  empirical  relation  shows  that  the  stiffness  and  damping  parameters, 
and  So  t  oan  be  established  with  the  me asurement  of  only  one  variable, 
which  proves  to  be  the  phase  angle,  rrtf. / £  „  If  aircraft  weight  and 

velocity  is  held  reasonably  oonstant,  it  becomes  exceptionally  convenient  to 
find  So  and,  JC0  from  the  following  rearrangement  of  the  equation  defining 
the  angle,  <L>  % 


d  Xjl.  (24) 

'  &o  /bo 


It  is  believed  that  the  normal  acceleration  response,  both  amplitude  and 
phase,  was  determined  with  good  accuracy,  and  therefore  it  is  possible  to  cheok 
the  previously  obtained  values  of  So  and  4  in  a  reliable  manner  without 
resort  to  the  pitching  velocity  measurements  in  any  way. 

In  the  actual  analysis  performed  in  reference  (3) ,  the  parameters,  So 
and  ,  were  calculated  for  each  frequency  test  point  in  order  to  allow  for 
the  variation  in  aircraft  weight  and  velocity.  Examination  of  the  results 
showed  that  the  scatter  in  the  values  of  S0  and  A?e  was  considerably  larger 
than  the  variation  in  speed  and  weight  encountered  in  flight.  The  experimental 
scatter  is  seen  to  be  a  function  of  the  accuracy  of  phase  angle  measurement. 

This  result  serves  to  emphasize  the  point  that  phase  angle  measurements  are  as 
important  or  even  more  important  than  amplitude  measurements.  The  resultant 
scatter  in  Sa  and  ST0  was  likewise  reflected  throughout  the  remainder  of  the 
analysis,  since  these  two  parameters  are  basic  to  the  calculation  of  the  rest  of 
the  derivatives.  Subsequent  computations  based  on  equations  (23)  and  (24) 
proved  that,  regardless  of  the  scatter  with  frequency,  the  originally  deter¬ 
mined  values  of  £>0  and  are  substantially  correct* 

During  this  study,  it  was  found  that  a  mechanical  harmonic  analyzer 
could  be  used  to  significantly  improve  She  accuracy  with  which  amplitude 
ratios  and  phase  angles  are  measured  on  the  oscillograph  records.  If  low  fre¬ 
quency  harmonic  distortion  is  present  in  the  reoord,  it  was  formerly  necessary  to 
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fair  a  sine  wave  through  the  trace  e  Fairing  was  similarly  required  when  a 
certain  amount  of  high  frequency  noise  was  encountered  in  the  oscillation  data. 
This  fairing  was  done  by  eye  and  necessarily  introduced  certain  errors  into  the 
data.  By  using  a  harmonic  analyser  to  select  the  fundamental  frequency  of  the 
input  and  output  trace ,  an  accurate  amplitude  ratio  and  phase  angle  can  be 
obtained  and  the  necessity  for  arbitrary  fairing  procedures  is  thereby  elim¬ 
inated. 

It  should  be  pointed  out  that  the  pitch  rate  gyro  used  in  the  F-80A 
longitudinal  flight  test  produced  an  output  trace  whose  wave  form  was  appre¬ 
ciably  distorted  with  additional  harmonics  and  at  tines  with  high  frequency  hash 
as  well.  These  performance  characteristics  made  it  very  difficult  to  obtain  a 
consistent  calibration  of  the  pitch  rate  channel  and  accordingly  the  accuracy  of 
the-  pitching  velocity  measurements  must  be  considered  somewhat  questionable. 
Indication  that  there  was  a  consistent  error  in  the  pitching  velocity  phase 
angle  measurement  is  given  by  figure  51  where  &&  is  shown  plotted  versus  the 
non-dimensional  frequency*  oof  0  From  equation  (12),  we  note  that 


f  -  -fc  <jn 


CO  f 


gcof  ^yc 


Hence* 


Since  and  4?  are  constants,  the  term  ooT-tasi&ttf should  be  constant 

with  frequency.  The  points  shown  on  figure  31  are  representative  of  the  results 
obtained  at  other  Mach  numbers  and  indicate  a  consistent  trend  which  is 
probably  caused  by  a  small  phase  error  in  the  dynamic  calibration  of  the  pitch 
rate  ohannel.  Numerical  studies  indicated  that  static  sensitivity  errors  in 
the  pitch  rate  measurement  would  produce  opposite  effects  on  the  extraction  of 
the  derivatives,  /7?£  and  .  For  example,  if  the  pitch  rate  were  assumed  to 
be  20  percent  greater  than  the  correct  theoretical  value,  an  approximate  16 
percent  decrease  in  and  24  percent  increase  in  ,  would  be  obtained  in 
the  reduction  to  derivative  form.  Note  that  the  F-80A  analysis  did  produce 
results  for  and  which  possess  this  qualitative  trend  with  respect  to 
the  wind  tunnel  evaluation  of  thess  parameters.  In  any  event,  this  investi¬ 
gation  was  somewhat  over-simplified  since  the  response  of  the  pitch  rate  channel 
was  not  perfectly  flat  over  the  frequency  range  covered  in  the  flight  tests. 

A  recapitulation  is  now  in  order.  The  flight  measured  value  of  static 
stability  is  in  drastic  disagreement  with  the  wind  tunnel  results,  with  the 
flight  value  being  approximately  one  hundred  percent  greater  at  0.3  Mach  number. 


WADC-TR-52-130 


-  28  - 


Sine©  the  original  F-80A  analysis  showed  the  tail-off  static  stability  to  be  in 
substantial  agreement  with  prediction,  the  increase  in  stability  was  explained 
by  a  stabiliser  effectiveness  which  was  considerably  larger  than  the  effective¬ 
ness  obtained  in  the  wind  tunnel.  Aerodynamic  theory,  however,  refuses  to 
acknowledge  the  possible  existence  of  these  unusually  large  values  of  and 

therefore  the  flight  measured  parameters  shown  on  figures  28  and  29  should  be 
regarded  as  extremely  questionable o  The  only  other  possible  explanation  for 
the  flight  observed  values  of  stiffness  is  that  (contrary  to  the  original 
analysis)  the  tail-off  static  stability  is  different  from  -chat  measured  in  the 
wind  tunnel.  Further  examination  of  the  problem  proved  that  this  was  the  ease. 

The  method  of  evaluating  the  tail-off  derivatives  using  the  responses 
to  a  tail  load  input,  was  outlined  in  the  previous  section.  In  the  original 
analysis  performed  in  reference  (S),  an  attempt  was  made  to  correot  for  the 
apparent  errors  in  the  tail  load  measurements,  but  in  view  of  the  results  that 
have  been  obtained  in  this  investigation,  it  is  believed  that  the  derived 
correction  procedure  was  inadequate  for  the  proposed  task.  No  attempt  was  made 
in  this  study  to  duplicate  the  original  techniques  developed  for  analyzing  the 
responses  to  tail  loads.  Instead  the  equation  of  motion  method  has  been  used, 
since  it  -gives  a  better  physical  insight  into  the  problem, 

Schumacher' s  curve -fi tting  technique^  (i,e,  equations  of  motion  method) 
was  first  applied  to  the  F-80A  response  data  to  determine  whether  agreement 
would  be  obtained  between  the  two  analysis  methods.  Computations  were  made  on 
the  basis  bf  three  variables  and  also  on  the  basis  of  four  variables,  in  the 
hope  that  the  tail-off  damping  in  pitch  oould  be  determined.  The  neoessary 
calculations  showed  that  could  not  be  extracted  from  the  response  data 

because  it  was  obtained  as  a  small  differenoe  of  two  large  quantities,  Values 
were  obtained  for  Mi  and  /??*,  ,  in  the  three  variable  solution,  which  agree 
exceedingly  well  with  the  results  obtained  previously.  These  computed  results 
are  shown  plotted  on  figures  28  and  29  for  0.3  and  0.7  Mach  number.  The  damping 
in  pitch,  ,  was  found  to  be  quite  large  in  this  three  variable  solution, 
but  in  view  of  the  assumption  that  ,  it  is  believed  that  some 

error  is  introduced  into  the  determination  of  this  quantity.  The  derivative,  /?&.  , 
is  also  plotted  on  figure  28.  .  * 

Note  that  this  equation  of  motion  analysis  oan  also  be  applied  to  the 
tail-off  equations  of  motion  as  given  by  equations  (13)  and  (14),  In  this  case, 
values  of  W*,  ,  rrfp.,  ,  and  could  be  determined  without,  the  necessity  of 

assuming  a  value  for  fCi/stoi  ,  since  Pol  is  not  a  variable  and  no  linear 
dependency  is  present.  This  analysis  scheme  was  first  tried  with  theoretical 
data  in  which  the  tail-off  damping  in  pitch  was  assumed  to  be  zero.  Again  it 
was  found  that  is  a  small  difference  term  and  even  when  using  supposedly 
perfect  data  a  value  of  /??.  -  -/  was  obtained.  This  result  re-emphasizes 

the  difficulty  with  which  Mp,  is  to  be  extracted  from  flight  data.  The  tail- 
off  static  stability,  however,  is  obtained  very  conveniently  by  this  method  and 
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of  course  will  possess  the  same  accuracy  as  -that  inherent  in  the  tail  measure¬ 
ments.  On  examination  of  equation  (14)  it  is  seen  tha't  for  a  given  total 
static  stability,  an  increase  in  tail  load  (meaning  an  increased  stability 
contribution  from  the  tail)  will  result  in  larger  values  of  negative  (unstable) 
static  stability  being  obtained  for  the  tail-off  airplane.  Since  the  analysis 
performed  in  reference  (?)  produced  an  extremely  large  value  for  the  slope  of 
the  tail  lift  curve s  measured  tail  loads  should  have  been  correspondingly 
higher  than  the  theoretical  predict! on.  Just  the  opposite  was  found,  however. 
It  therefore  becomes  necessary  to  show  that,  when  using  the  measured  tail  load 
data  or  a  calculated  load  based  on  the  actual  airplane  responses,  it  is  im¬ 
possible  to  justify  the  large  positive  values  of  ,  obtained  at  II  =  0.3  in 

the  original  analysis. 


In  performing  the  tail-off  equation  of  motion  analysis,  the  tail 
loads  computed  in  examining  discrepancies  between  predicted  and  measured  tail 
load  have  been  utilized.  These  tail  loads  were  based  on  the  measured  values  of 
( ^4.  ,  normal  acceleration  and  pitching  velocity.  It  was  necessary,  however,  to 
assume  values  for  ^  821(1  wil1  be  recalled  that 


wind  tunnel  values  were  used  for  the  first  two  parameters  and  a  flight  measured 
value  for  .  On  this  basis  the  horizontal  tail  load  was  calculated  for  all 
llaoh  numbers  at  the  20,000  feet,  27  percent  C.Go  flight  configuration.  Using 
the  tail  load  computed  for  0.3  Mach  number  (as  presented  in  figure  18) 
was  found  to  equal  0.411  or  ■  0.0853.  This  value  of  tail-off  static 

stability  is  shown  plotted  on  figure  28.  Note  that  the  difference  between  this 


point  and  the  estimated  tail-off  stability  is  approximately  equal  to  the  total 
increase  in  stability  exhibited  by  the  airplane  in  flight  over  that  indicated 
by  the  wind  tunnel.  At  a  Mach  number  of  0.7,  using  the  computed  tail  load  shown 
in  figure  20,  a  value  of  />?*/,  *  0.5061  or  \  ,  -  0.084  was  obtained. 

/ 1 

This  numerical  value  is  in  agreement  with  that  previously  determined,  but  again 
we  note  that  an  increase  in  stability  is  accompanied  by  a  similar  decrease  in 
the  instability  of  the  tail-off  configuration.  In  summary,  it  is  seen  that  the 
static  stability  contribution  of  the  tail  now  varies  with  Mach  number  in  a 


reasonable  manner.  Furthermore  ,  the  magnitude  of  ^ 


agreement  with  wind  tunnel  results. 


Cl, 


The  breakdown  of  inbo  values  of 


is  in  substantial 

and  is  the 


final  step  required  in  the  analysis.  Since 
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we  see  -that  the  large  value  of  obtained  in  the  original  analysis,  at 

0*7  Mach  number  for  example,  is  counteracted  by  a  large  value  of  »  It 

is  apparent  that  the  results  shown  in  figure  29  for  and  are  in 

error.  Given  values  of  the  transfer  function  constants,  and  ,  it  can 

be  demonstrated  that  /y?„  and  must  be  known  before  (Z+r  and  aZf 

can  be  correctly  evaluated.  On  solving  the  two  equations  (17)  and  (18),  it  is 
found  that 


Csr?.  - 


_ _ 


^7^^  *  \ 


CA 


The  above  expression  can  be  simplified  without  any  significant  loss  in  accuracy 
to  the  following* 


C<s  j&o  _  ,  /< 


£  47t/c 


On  inserting  typical  numerical  values,  we  find  that  the  two  terms  on  the 
extreme  right  are  very  important.  In  the  original  analysis,  /77^  was  generally 
too  large  (positive)  and  was  assumed  to  be  zero.  Both  factors  caused 

the  solution  of  this  equation  to  yield  values  for  (Tm £  which  were  too  large. 
In  the  equation  for  *  vis.  *  * 


j£-  = 


- 1  -t- 


it  can  likewise  be  shown  that  is  far  from  being  a  negligible  quantity, 

Values  for  can  also  be  obtained  from  a  solution  of  the 

below  two  expressions  for  and  ^  * 

-t 


WADC-TR-52-130 


-  31  - 


whore 


and 


=  Mp  ' 

=  m„  -m*, 


Again  it  is  Been  that  the  tail-off  derivatives  must  be  known  in  order 


to  accurately  compute  the  magnitude  of  and 


Unfortunately,  the 


extraction  of  ^p,  from  flight  data  does  not  appear  to  be  feasible  with 
presently  available  techniques.  No  serious  attempts  have  been  made  to  com¬ 
pletely  re-analyze  the  flight  data  presented  in  reference  (3),  since  the  required 
accurate  tail  load  data  are  not  available.  For  the  purpose  of  this  study,  it 
was  deemed  sufficient  to  establish  the  causes  of  the  subject  discrepancies  and 
to  indicate  methods  for  their  removal  rather  than  present  a  revised  numerical 
analysis  of  the  F-80  flight  data. 


4.  .  An  Equation  of  Motion  Analysis  Incorporating  a  Complex  Downwash. 

Statler  demonstrated,  in  reference  (l),  that  a  computed  response  will 
be  essentially  the  same  if  either  all  the  derivatives  are  treated  as  frequency 
dependent  or  just  the  downwash  is  assumed  to  be  frequency  dependent.  The  net 
result  is  such  -that  the  complex  downwash'  appears  to  cause  more  damping  than 
what  would  be  predicted  by  steady  flow  theory.  Computations  made  in  Appendix  C 
of  reference  (3)  indicates  that  the  frequency  dependent  damping  constant,  ^  , 
is  generally  higher  than  the  steady  flow  value  of  .  This  is  particularly 
true  at  higher  Mach  numbers.  Therefore,  if  tiier  e  is  any  evidence  to  show  that 
fee  complex  downwash  does  affect  the  response  in  accordance  with  the  theory 
present  in  reference  (l),  it  should  be  noted  as  increased  damping  over  and 
above  that  which  would  b9 .-predicted  by  the  simple  theory.  Values  of  , 


based  on  a  linear  analysis  of  the  F-80A  flight  data,  strongly  indicate  that 
this  evidence  exists.  Examination  of  figure  30  shows  that  the  measured  damping 
is  significantly  higher  than  prediction,  even  when  allowance  is  made  for  the 
tail-off  damping  which  was  neglected  in  the  original  calculations. 

On  the  basis  of  this  evidence,  let  us  assume  that  the  aircraft’s 
dynamics  are  more  accurately  represented  by  equations  of  motion  containing  the 
complex  form  of  downwash  rather  than  the  time  lag  oonoept.  What  are  the 
implications  of  analyzing  such  a  system  a&  if  it  were  non-complex?  Presumably, 
values  would  be  obtained  .for  />,  and  P'0  which  would  be  an  average  of  the 
actual  variation  with  frequency.  The  damping  parameter,  £>,  »  would  then  be 
larger  than  is  normally  defined  by  "steady  flow"  derivatives*  Accordingly,  we 
should  expect  that  the  experimentally  determined  lift  curve  slope  of  the  tail 
will  be  greater  than  that  indicated  by  the  wind  tunnel  for  steady  flow  con¬ 
ditions.  This  result,  of  course,  will  affect  the  experimental  determination  of 
all  the  basic  aerodynamic  parameters.  Based  on  the  data  which  have  been 
examined  in  this  study,  it  is  believed  that  the  above  situation  actually  exists 
and  has  influenced  to  an  uncertain  degree  the  analysis  of  the  F-80A  flight 
measured  responses. 

If  we  use  a  theory  which  consists  of  all  derivatives  constant  with 
frequency  except  /P*  (a  complex  )  and  which  is  devoid  of  the  so-called 

"time  lag"  term,  ,  it  appears  that  an  approximate  analysis  technique  may  be 
developed 'which  will  handle  the  introduced  oomplex  term.  In  accordance  with 
these  assumptions,  the  longitudinal  moment  equation  becomes 


D)oe  (25) 


where  f7 is  the  real  part  of 
and  is  the  imaginary  part  of 

*  // 

Figure  32  presents  a  plot  of  PI*,  and  P'1*/  as  a  function  of  the  reduced 
frequency,  °°  c  «  These  data  were  obtained  from  reference  (l).  If  the 
2  . 

-  /  // 

variation  in  P>*/  is  assumed  to  be  linear  with  frequency  and  the  curve  of  Ol*, 

is  assunmd  to  be  of  parabolic  form,  we  may  write 

mi,  =  A  +  &  IkVo 
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Substituting  -the  above  into  equation  (25)  we  obtain* 


f  Bco'7~)o{  +  (c.  +  £"  ^ 0)P G*  l 


Since  the  frequency,  c+ O  ,  appears  in  the  coefficients  of  the  variables,  it  ^ 
should  be  possible  to  obtain  a  unique  solution  for  the  unknowns  /),/?,  C  ,  £  , 
/7?a  ,  and  using  Schumacher’s  equation  of  motion  approach,,  Theoretically 

this  method  permits  an  approximate  determination  of  a  complex  by  means  of 
the  constants  A  »  £  ,  C  ,  and  £  0  Practically,  it  may  be  found  that  the 
constants  B  and  £  may  be  subject  to  the  "small  difference  difficulty"®  It 
is  seen  that  if  we  must  assume  B  and  equal  to  zero,  we  again  have  the 
equivalent  of  the  simple  system  with  a  time  lag  downwash  term,  where 

/9  -  mw 

and 


No  efforts  have  been  expended  to  apply  this  proposed  analysis  scheme, 
since  it  was  felt  that  the  available  time  should  be  spent  on  the  pitfalls  and 
weaknesses  of  the  analysis  of  a  non-complex  linear  system „  The  above  analysis 
method- is  presented,  primarily,  to  indicate  the  nature  of  the  simplifications 
required  to  transform  the  complex  /?9W  derivative  into  derivatives  which  are 
independent  of  frequency® 
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COMMENTS  ON  FUTURE  AEROELASTIC  RESEARCH 


The  subject  of  aeroelasticity  can  logically  be  divided  into  three 
separate  topics.  The  first  of  these  might  properly  be  called  "static"  aero- 
elastic  effects  and  would  include  such  items  as  loss  and  reversal  of  aileron 
control,  wing  and  wing-aileron  divergence,  loss  and  reversal  of  elevator 
control,  and  stabilizer-elevator  divergence,  and  the  less  spectacular  but 
sometimes  important  effects  of  deflection  on  various  derivatives.  These  phen¬ 
omena  may  be  observed  on  stiff  as  well  as  elastic  aircraft  and  are  primarily 
functions  of  the  dynamic  pressure  or  the  aerodynamic  loading.  The  second  topic 
is  the  well  known  phenomenon  called  "flutter"  which  results  from  the  dynamic 
coupling  existing  between  the  unsteady  aerodynamic  farces  and  the  combined 
structural  bending  and  torsion  modes  of  an  aerodynamic  surface  or  control.  This 
phenomenon  is  a  function  of  true  airspeed  and  the  natural  frequency  of  the 
structural  flutter  mode  and  generally  occurs  at  frequencies  far  above  the  natural 
frequency  of  the  aircraft's  rigid  body  motions.  The  third  topic  consists  of 
the  relatively  newer  phenomenon  in  which  the  structural  dynamics  have  been 
observed  to  significantly  affect  the  rigid  body  response  of  the  airplane.  This 
aeroelastic  effect  is  a  function  of  the  natural  frequencies  of  the  aircraft, 
structure  and  increases  in  magnitude  as  these  frequencies  approach  those  of  the 
aircraft  rigid  body  response. 

Although  the  influence  of  aeroelasticity  on  the  dynamics  of  the  F-80A 
was  quantitatively  investigated,  it  should  be  realized  that  this  configuration 
does  not  provide  an  example  of  the  third  named  aeroelastic  effect  in  which 
structural  motions  are  found  to  couple  dynamically  with  the  motion  of  the  air¬ 
plane's  center  of  gravity.  Since  the  F-80A  is  a  structurally  stiff  airplane, 
it  is  seen  that  aeroelastic  effects  are  essentially  scalar  rather  than  dynamic* 
Note,  also,  that  the  structual  motion  (fuselage  bending)  is  very  small  as  is 
indicated  by  figure  26.  These  considerations  lend  weight  to  the  argument  that 
the  F-80  airplane  is  not  a  suitable  aeroelastic  research  test  vehicle. 

As  aircraft  continue  to  become  larger  and  more  flexible,  coupling 
between  the  structural  and  rigid  body  modes  will  naturally  become  more  pronounced. 
This  trend  has  already  been  dramatically  demonstrated  by  difficulties  en¬ 
countered  in  present  day  aircraft  in  which  structural  modes  have  coupled  with 
various  types  of  automatic  control  systems.  The  problem  th*n  arises  of  being 
able  to  theoretically  predict  the  response  of  the  flexible  airplane  which  has  a 
large  number  of  structural  degrees  of  freedom.  This  situation  is  further  com¬ 
plicated  by  the  introduction  of  the  swept  wing  since  wing  bending  then  produces 
changes  in. the  span-wise  load  distribution.  These  facts  would  indicate  that 
the  demonstration  of  the  feasibility  of  accurately  predicting  the  response  of 
to  "aeroelastic"  airplane  is  certainly  in  order.  Such  a  demonstration  should 
certainly  involve  the  full  scale  flight  tests  of  a  truly  aeroelastic  aircraft. 
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COHCLUSIOSS  ASP  SECOMMBHDAT IOHS 


On  completion  of  the  study  described  herein,  it  has  been  concluded 

thats 

1.  Whereas  the  determination  of  the  short  period  longitudinal  stiffness 

and  damping  parameters  from  dynamic  flight  data  is  a  straightforward  task,  all 
of  the  basic  aerodynamic  quantities  cannot  be  determined  (including  etc /ctoi- 
and  )  tinless  means  are  available  for  determining  the  tail-off  static 

stability  and  damping  in  pitch. 

2,  The  tail  load  discrepanoy  observed  in  reference  (5)  is,  to  a  large 
degree,  due  to  an  error  in  measurement  due  to  limitations  in  possible  strain 
gage  locations. 

5.  the  effects  of  fuselage  bending  and  elevator  twist  an  the  response 
of  ihe  F-80A  are  predominantly  scalar  and  therefore  difficult  to  perceive  in 
the  experimental  measurements. 

4.  The  lack  of  agreement,  shown  in  reference  (3),  between  predicted  end 
measured  responses  at  low  Uaoh  numbers  is  due  to  prediction  being  based  on  a 
wind  tunnel  evaluation  of  tail-off  static  stability  which  is  different  from 
that  obtained  in  flight. 

6.  From  the  response  prediction  standpoint,  the  errors  introduced  by 
the  estimation  of  the  first  order  aerodynamic  derivatives  are  large  enough  to 
completely  mask  the  effeots  of  the  seoond  order  terms,  the  unsteady  flow  effects, 
sad  the  possible  aeroelastio  effects  for  the  F-80A  airplane. 

6.  The  numerical  results  obtained  in  the  original  F-80A  analysis  for 
stabiliser  effectiveness  and  slope  of  the  dowswash  are  in  error  due  to  inoorreot 
values  of  tail-off  static  stability  and  tail-off  damping  in  pitch  being  used 
in  the  analysis. 

70  There  is  evidenoe  that  the  measured  damping  of  the  F-80A  airplane  is 
signifioantly  higher  than  what  would  be  predicted  by  linear  theory,  giving 
indication  that  the  complex  dow awash  does  effeot  the  response  in  accordance 
with  the  calculations  made  in  reference  (l). 

8.  The  most  significant  dynamic  aeroelastio  problem  at  this  is  the 
increasing  ooupling  between  structural, automatic  control  system,  and  rigid 
body  modes  resulting  from  present  day  aircraft  design  trends  toward  more  flex¬ 
ible  aircraft. 

9.  The  F-80A  airplane  is  not  a  suitable  airplane  for  examining  dynamic 
aeroelastio  effects. 

Additional  or  secondary  conclusions  may  be  stated  as  follows; 
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lOo  Since  careful  estimation  of  the  static  stability  is  necessary  for 
aocurate  prediction  of  the  longitudinal  response,  it  would  seem  advisable  that 
static  flight  tests  be  made  to  verify  wind  tunnel  results. 

11.  In  the  analysis  of  tail-on  longitudinal  response  data,  the  damping 
and  spring  constants,  /},  and  ,  are  defined  by  the  measured  phase  angle, 

P  *>}/£  .  A  statio  sensitivity  error  in  the  elevator,  pitohing  veloeity,  and 
normal  acceleration  traces,  therefore,  affects  only  the  extraction  of  the 
derivatives,  and  yWg  . 

12 o  A  mechanical  harmonic  analyzer  can  be  used  to  significantly  improve 
the  accuracy  with  which  amplitude  ratios  and  phase  angles  can  be  measured  from 
sinusoidal  oscillation  data  in  which  a  disturbing  amount  of  harmonic  distortion 
is  present. 

13.  It  does  not  appear  possible  to  extract  the  tail-off  damping  in  pitch, 

,  from  longitudinal  response  data  with  any  of  the  analysis  techniques 
tried  to  date. 


The  investigation  reported  on  herein,  together  with  the  above-stated 
conclusions,  has  led  to  the  following  recommendations. 

1.  Further  studies  should  be  made  to  develop  methods  which  would  simplify 
and  expedite  the  measurement  of  horizontal  tail  loads.  These  studies  should  be 
accompanied  by  the  development  of  additional  analytical  techniques  for  deter¬ 
mining  the  tail-off  stability  derivatives. 

2.  Efforts  should  also  be  directed  toward  analyzing  longitudinal  response 

data  on  the  basis  of  equations  of  motion  consisting  of  all  derivatives  constant 
with  frequency  except  (a  complex  *££ /f^oL  )  and  which  are  devoid  of  the 

so-called  "time  lag"  term,  • 

3.  The  feasibility  of  accurately  predicting  the  response  of  an  aero- 
elastic  airplane  should  be  demonstrated  by  means  of  full  scale  flight  tests 
with  a  truly  aeroelastic  aircraft. 
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